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I. INTRODUCTION 

The two dimensional form of sp 2 hybridized carbon, 
graphene, is a flexible one-atom-thick soft membrane em- 
bedded in three dimensional space. Due to an anhar- 
monic (nonlinear) coupling between bending and stretch- 
ing phonon modes^, it comes with ripples whose typical 
height scales with sample size L as L* , with t = 0.6. 
Besides providing structural stability, ripples also play 
an important role in graphene's electronic properties^. 
Graphene has two atoms per unit cell, which results in 
two conical points K and K' per Brillouin zone where 
band crossing occurs^, see Fig. [TJ Near these crossing 
points, the electron energy is linearly dependent on the 
wave vector, see Fig. [2j This behavior is robust with 
respect to long-range hopping processes beyond simple 
nearest-neighbor, tight-binding approximation^ because 
it follows from symmetry considerations^. Here we con- 
sider the constrained relativistic problem for the carriers 
in two dimensions and discuss the effects that arise from 
the emergence of a geometry induced potential Vq on the 
electronic properties. 

The effect of random surface curvature ripples on the 
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FIG. 1: Left: The lattice structure of graphene is comprised 
of two interpenetrating triangular lattices A and B with lat- 
tice unit vectors ai and a2. Right: The Brillouin zone where 
the Dirac cones are located, K and K'. The reciprocal lattice 
vectors are bi and b2. 



electronic properties of monolayer graphene is primar- 
ily described via gauge potentials^. However, the effect 
of geometry induced potential a la da CostcP in mono- 
layer graphene has not been investigated. The effect of 
geometric potential for bilayer graphene is discussed in 
Ref. [5] but in this case the quantum problem is not rela- 
tivistic. Indeed, a rigorous derivation of the effect curva- 
ture ripples have on the spectrum of the surface carriers, 
electrons and holes, confined on monolayer graphene re- 
quires a consistent confining procedure for massless Dirac 
fermions onto a surface embedded in three dimensional 
space, which we present here. 

Low dimensional systems are realized by constraints 
(holonomic or nonholonomic) that classically reduce the 
degrees of freedom available for the particle, namely in 
two-dimensional electron (or hole) systems (2DES), typ- 
ically realized in III-V semiconductor hetero junctions 9 . 
The barrier potential at the heteroj unction interface, cre- 
ated by differing electron affinities, leads to a clear sepa- 
ration of energy scales for the motion in the 2D plane and 
motion along the direction normal to the plane. Thus the 
particle's wave function is separable into normal and sur- 
face components and, at low energies, dominated by the 
surface component!^. Constraints can thus be realized 
by introducing external potentials forcing the system to 
occupy less degrees of freedom^]. 

The quantum properties of a nonrelativistic parti- 
cle constrained to an arbitrary orientable surface are 
knowrPl the geometry of the surface induces an attrac- 
tive geometric potential in the Schrodinger equation in 
curvilinear surface coordinates 

h 2 

VdaCosta (91,^2) = -7— («1 - , (1) 
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where m is the effective mass of the particle, h is Planck's 
constant, (91,(72) denote surface coordinates, and k 1; k 2 
are the two position-dependent principal curvatures of 
the surface. This potential is purely a result of particle 
confinement, and is independent of the electric charge of 
the particle; it is therefore the same for electrons and 
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holes. This result is applicable in the limit 
e K — > 0, 



(2) 



where eo is the thickness of the two-dimensional surface 
and 



k = max{fti, K2}. 



(3) 



is the surface curvature. Here eo will correspond to the 
width of the normal to the surface quantum well in 2DES 
where particles are confined. However, absent a truly two 
dimensional system that can be easily bent, the effect of 
geometric potential on the electronic band structure has 
been justifiably ignored in device engineering up to now. 
Graphene represents a unique material which can exhibit 
the effects produced by the geometric potential. 
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2: Left: The electronic structure of graphene^ 
: ±iy / 3 + /(fc) - t'f(k), where f(k) = 2 cos (V3k y a) + 



4cos (^^k y a j cos (|fc^a). Here k = (k x ,k y ), t 



= 2.8 e\^, 

t' = 0.1 eV^ and a = 1.42 A. The plus sign applies to the 
upper (tt) and the minus sign to the lower (71*) band. Right: 
The band structure zoomed in on one of the Dirac points 

where E± = ±hvF jfcj . The Fermi velocity vf is constant and 
the electron-hole symmetry is not broken if next-to-nearest 
neighbor hopping is neglected. 



Curving graphene can have three major microscopic 
effects: (i) changing the distance between carbon atoms, 
which modifies the hopping amplitude between neighbor- 
ing sites and is, in general, very costly due to the large 
spring constant of graphene w 57 eV/A^, (ii) a rota- 
tion of the p z orbitals and (iii) rehybridization between 
7r and a orbitals 13 . Curvature in sp 2 hybridized carbon 
allotropes also modifies the Raman spectra through ex- 
erting strain in the underlying latticed The rotation 
between the orbi tals ca n be understood within the Slater- 
Koster formalis m 1 13 ! 15 !. When tt orbitals are not parallel 
the hybridization between them depends on their rela- 
tive orientation, which is a function of local curvature. 
Furthermore, rotation leads to rehybridization between 
tt and a orbitals which shifts the energy even further. 



As a result, Dirac fermions are scattered by ripples and 
curved regions through a geometric potential Vq which is 
derived in Section [Hi] Section [TT] discusses the geometric 
mechanism which shifts the Fermi energy thus leading 
to the emergence of p-n junctions in monolayer grap hene 



IV deals 



akin to the case of bilayer graphene 21 . Section 
with periodic corrugations which produce an energy band 
spectrum. In Section[V]we pose the inverse problem, that 
is how to engineer surfaces with prescribed quantum be- 
havior. In Section [VI| the membrane aspect of graphene 
is discussed in view of the constrained two dimensional 
Dirac equation derived in Section |III[ In Section VII we 
summarize our main findings. 



II. p-n JUNCTIONS 

In this paper, we argue that due to its essentially 
two-dimensional nature e n 0, gapless linear band 
structure, and exceptional material strength, monolayer 
graphene presents an excellent candidate for exploration 
of two-dimensional systems with curvature. Dimensional 
analysis implies that for graphene^, the geometric poten- 
tial Vq must scale as 
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(4) 



where vf ~ c/300 is the Fermi velocity of massless car- 
riers in graphenelSI. Here c is the velocity of light in 
vacuum. 

Now, let us discuss the applicability of a constrain- 
ing procedure to monolayer graphene. To map a curved 
monolayer graphene sheet on a relativistic quantum prob- 
lem for a particle on a curved surface, two criteria must 
be satisfied^ 



and 



/cfog -C 1 



kclc -C 1, 



(5) 



(6) 



where ag = 1.42A is the carbon-carbon distance in 
graphene, \zf is the Fermi wave vector, and k is the sur- 
face curvature given in ([3]). The first criterion ensures 
that monolayer graphene can be treated in the continuum 
limit and the detailed lattice structure can be ignored. 
For monolayer graphene near the Dirac point where the 
relativistic description holds hp — and the first crite- 
rion is always satisfied. The second criterion ensures that 
the detailed lattice structure is not perturbed on the cur- 
vature length-scale and the relativistic description holds 
true because the hexagonal symmetry is not disrupted; 
and because the thickness of monolayer graphene sheet 
is e ~ 3ac it also ensures that e n -C 1, or that the 
monolayer graphene sheet with nonzero thickness can be 
considered a surface. Indeed, for carbon nanotubes with 
typical radii rj 10 — 30 A and for typical surface rip- 
ples 1/k w 100 A on monolayer graphene^, the second 
criterion is satisfied too. 
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As a result, for monolayer graphene k ^> kp = the 
geometric potential Vq becomes dominant over the Fermi 
energy Ep. For example, let us consider graphene with 
n-type carriers in the flat region where Ep = 0, then 
the geometric potential will lower the Fermi energy from 
the conduction band to the valence band, Ep — — \Vq\- 
Note that as long as the geometric potential is constant 
or slowly varying, as is the case for ripples, this result 
is robust and is independent of any quantum numbers 
the carriers may have (chirality and/or electric charge). 
Therefore, the curved regions will have p-type carriers 
with a natural p-n junction created between the curved 
and flat regions. For the surface z = f(x), the geometric 
potential will create p- and n-type strips where the loca- 
tion, width, and polarity of each strip is determined by 
the local curvature dependent Fermi energy Ep(x), see 
Fig. [3] and the derivation which follows. 
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FIG. 3: n-p-n junction (transistor) formed by bending 
monolayer graphene shown in cross-section. The Fermi en- 
ergy is a function of the position and therefore of curvature 
Ef{x) = — I Vg I oc —hvF-ji along the sheet. 



III. CONSTRAINING THE DIRAC EQUATION 

Let us consider a corrugated graphene sheet along one 
of the coordinates, for example x. We will demonstrate 
the confining procedure leading to the appearance of a 
curvature dependent potential in the Dirac equation de- 
scribing the carriers, electrons and holes, in graphene. 
The confining procedure starts with the introduction of 
an external to the graphene sheet confining potential 
Vj_ which forces the quantum system to remain on the 
two dimensional surface thus simulating a mechanical 
constraint. This potential for graphene stems from the 
Coulomb interaction since as carriers leave the surface of 
the graphene sheet (their wavefunction is not strictly two 
dimensional) they distort the charge balance and a po- 
tential arises which forces them back on the surface. This 
potential is similar to the potential with which a charged 
plate acts on a single charge above it. The properties of 
V± will be introduced later. 



The curved two dimensional graphene surface em- 
bedded in a flat euclidean three dimensional space is 
parametrized as follows 



f{x, y) = xe x + ye y + f(x)e z , 



(7) 



where (e x , e y ,e z ) is the cartesian coordinate system asso- 
ciated with the embedding space, that is the laboratory 
frame. 

In order to derive a moving coordinate system 
(ei,e*2,e3) associated with the surface at each point we 
need to differentiate the parametric equation Q 



ei 



e-2 



e3 



dr 
dx 

dr 

dy 
dr 

dx 
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1 dr = 

dx ~ ^1 + f- 



r 



dr dr 
dx dy 



, (8) 
(9) 
(10) 
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The vector e*3 is the normal to the surface at each point. 
The other two vectors e\ and ei span the tangent space at 
each point. Hereafter /' = df /dx denotes differentiation 
with respect to x. The moving frame (e\, 62,63) can be 
related to the immobile laboratory one 
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where 
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From the above two formulas we establish two important 
relations to be used later in the calculations 
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Next we introduce a coordinate system of the embedding 
space associated with the corrugated graphene sheet Q 



R(x,y,e) = r(x,y) + ee 3 (x). 



(15) 



Here e measures the deviation from the surface in normal 
direction along e*3. The metric tensor of the embedding 
space can be computed according to 



G, 



dR dR 

di dj 1 



(16) 
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where 



dR 
dx 
dR 
dy 



= y/i + /' 2 ei + 



e2, 



dR 



dx 
e*3- 



A straightforward computation yields 



dx 



-/te 3 , 



de 3 
dx 



nei, 



(17) 
(18) 

(19) 



where hereafter the nonvanishing curvature k of the 
graphene sheet is 



f" 

K 1 + /' 2 ' 

and should not be mistaken for pi). 

In the mobile basis (e\, e*2, S3) the metric tensor is 



(20) 



G 




where 



A = v 7 ! + /' 2 + ■ 



(21) 



(22) 



In the moving orthonormal frame the gradient operator 
becomes 



V = e 1 
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<~-2 



d_ 
dy 



e 3 



d 

de' 



(23) 



With units h — vp — 1 the hamiltonian of the massless 
Dirac equation Hty = describing carriers, belonging 
to different sublattices (via the two component spinor 
wavef unction in corrugated graphene is H — — zdf.V. 
Here a is a vector consisting of 2 x 2 Pauli matrices. 



In the coordinate system (15) the Dirac hamiltonian 
takes the form 

_H _ _ 

H = -ia.ei — - ia.e 2 ^ ia.e 3 — . (24) 

A dx dy de 

The wave function \1/ is normalized according to 

Adxdyde = 1, (25) 



where t denotes hermitian conjugation. The factor A 
in the integration measure complicates calculations and 
obscures hermiticity. We define a new wavefunction 



iP = A 1 / 2 * 
which has the normalization condition 

tp'ip dxdyde = 1. 



(26) 



(27) 



The original Dirac equation can be recast in terms of 
an energy operator H — A 1 / 2 'HA _1 / 2 which acts on the 
wave functions %b 



H= —lot.e-i^ ia.e 3 — — \-i—-a.e 3 

dy 



de ' ~2A' 
d 1 0A N 
dx ~ 2Al)x 



(28) 



This hamiltonian is difficult to compute with because the 
matrices a.eTj are ^-dependent. The momentum opera- 
tors here are expressed with respect to the moving coor- 
dinate system Sj while the pseudo-spin operators a are 
expressed with respect to some fixed rectilinear frame. 
To remedy this we change the pseudo-spin basis through 
a unitary transformation f2(x) 



fl(x)a.e j (x)Q.^(x) 



(29) 



holding true for all x. Here the x-dependence of the unit 
vector Sj is made explicit. 

To begin with, the moving frame undergoes an evolu- 
tion as traversing the x direction 



dej_ 
dx 



for j = 1,2,3. 



Here 



e 3 x 



de 3 
dx 



-ne 2 



(30) 



(31) 



is a smoothly varying instantaneous "angular velocity" 
associated with the evolution of the frame. 

Next, we suppose Sl(x) — "Pexp | J x dx'"E.uj(x') 

where V denotes path ordering along x and E is to be 
determined. This yields the following equation for 



dfl 

dx 



-E.cJfi. 



(32) 



After differentiating (29) with respect to x, we arrive at 
dfi f 



dfl 

dx 



cJ.ejO' + Via. 



dx 



dej 
dx 



fi f =0. 



Now using ( 30 1 and ( 32 ) we obtain 

de 



E.cJ, a, 



dx 



(33) 



(34) 



The form we have chosen for f2 points to a change of 
pseudo-spin basis to spinors representing carriers in sub- 
lattice A and B (Fig. 1) in graphene's corrugated surface 
described locally by e\ and e 2 . 

We write down a small computation table 



E.w, ai = — 2ina 3 , 
E.c3, a 2 = 0, 
E.w, «3 = 2ina\. 
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From the middle equation follows E.cJ = Ca 2 - If we 
choose C = k the above table transforms into 

[a 2 ,ai] = -2ia 3 , 
[a 2 ,a 2 ] = 0, 
[a 2 >Q!3] = 2i«i. 

These relations coincide with the algebra of the Pauli 
matrices [a a , a b ] = 2ie abc a c a x a 2 = ia 3 . 

We define a new spinor wavefunction to take advantage 
of this transformation 



(35) 



in addition to a new Dirac hamiltonian H = f2_fff^ act- 
ing on x- Although this transformation replaces the a- 
matrices as intended, it generates a gauge term in H: 



i a ~~drt 1 - 
— -ila.ei— — = — — — (XiIj.uj. 
A dx 2A 

Using the above given tables we find 

- _ i _ _ dW in 



(36) 



a 3 . 



This term cancels the third term in ( 28 1 . Finally, the 



three dimensional Dirac hamiltonian is simplified 



H = 



8 



8 



= —iu 2 IQ.3 tt- 



'8y 



A 



8e 

d_ _ J_9A 
dx ~ 2A ~8x 



(37) 



The transformations we have performed reduce the 
original Dirac problem to the study of this perturbed 
Dirac hamiltonian in the space —eo/2 < e < eo/2, where 
eo ~> 0, between two parallel identical corrugated sur- 
faces < x < L x , < y < L y with respect to a fixed 
rectilinear frame (see Fig. [4]). This simplification gave 
rise to a variety of curvature dependent terms in ( 37 1 . 




FIG. 4: The gradual confining procedure eg — > and the 
domain where the Dirichlet boundary conditions apply for 
(371. 



Despite its two dimensional nature, graphene has three 
acoustic phonon modest The two in-plane modes (LA, 



TA) have a linear dispersion relation, whereas the out 
of plane mode (ZA) has a quadratic dispersion relation 
pointing to the possibility of having an effective mass 
term aim± included in the above hamiltonian which is 
three dimensional and contains information for the off- 
surface processes as well. As soon as the carriers leave 
the two dimensional world of the graphene sheet where 
symmetry precludes the mass, they attain mass. There- 
fore we will treat the transverse dynamics in the non- 
relativistic limit. Furthermore, according to the ideas 
presented in the beginning of this section, we shall now 
consider the idealistic spatial potential Vj_(Ao,e), where 
Ao measures the squeezing strength of this potential 



lim Vj_(A ,e) 



0, 



0, 



oo, e ^ 0. 



(38) 



For example, we can imagine a realistic harmonic binding 
V±(Xo,e) = ^m±X^e 2 with Ao going to infinity which 
yields (e 2 ) ~ h/m±\Q in the nonrelativistic limit. 
Next we look for eigenstates of (37) in the form 



x(x, y, e) = ^2 y) h *( e ) 



(39) 



where h\ (e) are eigenstates of the transverse one dimen- 
sional hamiltonian 

d 

H± = -ia 3 — + a im± + V ± (e). (40) 

The label A denotes the quantum numbers necessary 
to specify the eigenstates h\(e) of H± completely and 
physically speaking it is the "transverse energy" due to 
confinement. The boundary conditions are vanishing 
(Dirichlet type) at e = ±e /2. We need to take the non- 
relativistic limit when solving H± in order to overcome 
the well known jittery behavior (due to Heisenberg prin- 
ciple) of the wavefunction when confined in too small an 
interval ej^- 

The remaining two dimensional part amounts to 



H' 



-lOt-2 



d_ 
9y 



-ax 



d_ 

dx 



1 <9A 
2A&T 



(41) 



Since we have already factored out the transverse part, 
we take the limit e — > in (41 ). This action is in agree- 



ment with the standard procedure in constrained quan- 
tum dynamics in the case of small eo ~ 3acpl The small- 
ness of eo in graphene is combined with the relatively 
small curvature of the naturally occurring corrugations. 
Thus for graphene the following applies 



8k 
dx 



M < l, 

rendering the approximation 



< i, 



(42) 



(43) 
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justified. Indeed, 
. d 



H' 



-10L2 



9y 



I0i\ 



d , K (x) f 



Introducing a new variable, which is the line length along 
the corrugated surface 



Vl + f' 2 dx', 



instead of x and 



K[s(x)} = k(x) 



r 



we may rewrite H' as 



where 



H' = —ia-i— ia\ 

dy 



V G (s) 



-V G (s] 



-K(s). 



(44) 



(45) 



(46) 



(47) 



As is often the standard approach with the Dirac equa- 
tion, it is easier to work with the square of the Dirac oper- 
ator instead of the first order form H' 2 if}\ = (E—\) 2 ip\ — 
E^tpx, where E is the total energy and E e i is the energy 
of the electrons on the surface. We have 

d 2 



H' 2 = - 



ds 2 



IdK l R d 

2lh ~ 2 ds 



1 



dy 2 ds 2 2 ds 2~ ds 4 
Looking for solutions obeying the ansatz 



K 2 (s). 



(48) 



we simplify the squared hamiltonian acting on rj\ consid- 
erably 



H' 2 = - 



where 



V(s) 



d 2 ^ 
dy 2 



1 
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d 2 
ds 2 



dK 



V(s), 



4^+3K z {s) 
ds 



(49) 



(50) 



The hamiltonian ( 49 ) with a non- vanishing potential 



( 50 1 possesses the usual phenomenology of bound and 



scattering states. Moreover, band structure is allowed 
for periodic V(s). Remembering the geometrical origin 
of (50), we may state that all of the above mentioned 



instances for graphene are of geometrical origin including 
decoherence as a function of scattering and interaction 
through Vq. 



FIG. 5: A periodically corrugated surface f(x) — sin(a:) 



IV. ENERGY BANDS 

Let us explore a periodically corrugated graphene sheet 
(see Fig. [5] ) parametrized as follows 



f(x) = esin . 



(51) 



Suppose ie/ a) 2 <C 1 is a small parameter then the fol- 
lowing approximation holds true up to second order in 
this small parameter 1 + /' 2 « 1 and s ss x. 
For the potential V(s) we obtain 



V(x) = 



le 2 

~4a 4 



cos 



2x 



2 7 a 6 



le 2 (2x 

- - -r cos 
4 a 4 



Introducing the dimensionless variable z — s/a, the s- 
componcnt of the squared Dirac hamiltonian ( |49[ ) acting 
on a factorized wave function rj\(y,s) = 8(y)<fr(s) takes 
the form of the Mathieu equation 



dz 2 



a — 2q cos (2z) 



(z)=0, 



(52) 



i 

where a = ak s and q = — | ^ ■ Here k s is the wave 
number along the s-degree of freedom. The solutions 
of the above are given in terms of Mathieu functions^ 
and a band gap structure naturally emerges as a conse- 
quence of the properties of the Mathieu equation^. The 
y-component of the solution obeys the standard harmonic 
oscillator equation and will not be discussed here. 

According to the Floquet's theorem, the general solu- 
tion to the Mathieu equation is 



z) = cie iz h(z) + c 2 e~ liZ h(-z) 



(53) 



where h(z) is a periodic complex valued function with 
period 2tt and the characteristic exponent £ is a definite 
complex valued function of a and q. Here the constants 
Ci are arbitrary up to Cic 2 7^ 0. The characteristic ex- 
ponent £ = ip is imaginary for spatially undamped so- 
lutions, while it is real or complex for spatially damped 
solutionflH. Spatially undamped solutions exist in spe- 
cific regions only defined by £ = £[a(a, q s ), q(a, e)} and are 
refered to as allowed energy bands. There are four types 
of periodic solutions, even z-parity ce p (q, z) and odd z- 
parity se p (q, z) in two classes, even-integer at p = 2n 
and odd-integer at p = 2n + 1, with respective eigen- 
values (respectively eigenenergies) a p (q) to ce p (g, z) and 
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b p (q) to se p (q,z). The functions are normalized so as 



ce (q, z) w 1, 



1. At small \q\ 



1st 

8 a 2 



ce p {q, z) w cos (pz), se p (q, z) « sin (pz) 



« 1, 



(54) 



and a p ~ p 2 , 6 p ~ p 2 . These integer class periodic func- 
tions define the edges of the allowed energy bands. 

Mathieu eigenvalue equation defines allowed energy 
bands at q < (upper-sign branc h) and at q > 
(lower-sign branch) respectivel y^ 2 * 23 *. As q is increased, 
the allowed bands get narrower turning to single levels 
which is not the case here since q < 0. The allowed 
band edges in the energy versus q plane are described by 
Mathieu eigenvalues a p (q) and b p (q), or alternatively, by 
the corresponding periodic eigenfunctions ce p (q, z) and 
se p (q, z) at integer p (p = 1, 2, 3, . . .), while the solution 
for each energy in the interior of the band is generated 
by pairs of Mathieu functions at intermediate non- integer 
p. Whithin an allowed band, Mathieu eigenvalues and 
eigenstates are continuous functions of the wave number 
k = p7r, where < p < 1 for the first Brillouin zone^. 
At q < 0, which is the case here, the allowed bands are 
defined as follows: (ao, a%), (b±, 62), {0*2, 0-3), etc. 



V. INVERSE PROBLEM 

We showed that the relativistic quantum dynamics of 
massless carriers in corrugated graphene can be mapped 
on a two dimensional Schrodinger problem ( 49 ) with the 
geometric potential ( 50 1 . Now we pose the inverse prob- 



lem (for the nonrelativistic case it is discussed in Ref. 
|24j ) . namely we may impose a prescribed potential U(s) 
and solve for the graphene surface profile f[s(x)], i.e. we 
solve the following third order nonlinear ordinary differ- 
ential equation in terms of f(s) given some boundary 
conditions 



U[s(x)} = V[s(x)}. 



(55) 



However, here we will focus on a simpler problem. For 
the constrained two dimensional Dirac equation, where 
the geometric potential is Vg[s(x)], we can pose a similar 
inverse problem for a prescribed potential R(s) 



R(s) = V G (s). 



(56) 



Let us answer the question which surfaces do not produce 
a geometric potential, that is 



R( 8 )=0=>V G [8(x)]=0. 



Using ( 47 1 we can show that this is equivalent to 



(57) 



(58) 



The solution for the graphene profiles which do not pro- 
duce geometric potential corresponds to flat surfaces only 



f{x) = C\x + C 2 



(59) 



where C\ and Ci are arbitrary real constants. Any cor- 
rugation will generate a geometric potential Vg =/= 0. 



VI. AN ELECTRONIC MEMBRANE 

Graphene is an example of an electronic membrane apt 
to gating and probing 13 which adds a new perspective to 
the physics of membranes^. In this section we estab- 
lish a connection between the membrane aspect and the 
electronic properties of graphene through the geometry 
induced potential Vg and the Dirac equation constrained 
in two dimensions ( 46 ) . To build an effective theory we 
introduce the bending free energy Eb enc [ describing an 
almost flat membrane in the thermodynamic limit 26 



dq\dq 2 ^/g 



7 + ^(2MY 



vK, 



(60) 



where (91,92) are the surface's coordinates, g is the de- 
terminant of the metric, M and K, are the Mean and the 
Gaussian curvature, respectively. Here 7 is the tension, 
/i is the bending rigidity and v is the Gaussian rigidity. 
Besides /1 ~ 1 eW-ZI, 7 and v are presently not known. 
If we neglect the Van der Waals interaction between the 
graphene membrane and an external support, the result- 
ing free energy using formulas from Refs. |13l28j is 



E, 



bend 



J dxdy^l + /'■ 



1 2 

7+ 



(61) 



where k is given by ( 20 ) . The Gaussian total curva- 



ture term vanishes because it is a total derivative and 
for a membrane with a fixed topology it gives a con- 
stant (which is neglected) according to the Gauss-Bonet 
theorerrPH. Now, the above expression has to be com- 
bined with the following energy functional stemming 
from J39l and d46l) 



E, 



E 



dxdyde\hx(e)\ 2 tp{ 
1 d 



-id2 



d_ 
dy 



(62) 
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The coupled problem 8(Eb en d + E e i) — has to be solved 
self-consistently in order to give graphene's equilibrium 
corrugations, which will be a function of the confined 
eigenstates h\(e) and the strength of the confining po- 
tential V±(e). 



VII. CONCLUSION 

In conclusion, we have demonstrated that for n-doped 
graphene the geometric potential Vg due to confinement 
in the Dirac equation for the massless carriers dominates 
over the Fermi energy Ep in the bent regions (where the 
curvature is nonvanishing) and leads to the formation of 



alternating p- and n- type regions, that is p-n junctions 
of pure geometrical origin. This may have a practical 
application in constructing wavelength-specific (because 
Vg can be tuned) solar cells based on graphene. Next, 
we have derived proper constraining procedure leading 
to the two dimensional Dirac equation for the carriers 
in corrugated graphene. The geometric potential due to 
confinement Vg emerges as a consequence of the quan- 
tum mechanical property that the wavefunction is three 
dimensional and when constrained "experiences" the cur- 
vature of the sheet. The inverse problem, namely the 
engineering of graphene surfaces with prescribed quan- 
tum properties, is posed. It opens an avenue for device 
engineering based on corrugated graphene. Periodic cor- 
rugations lead to the emergence of band gap structure 
in the energy spectrum, which can be tuned by vary- 
ing the amplitude and the period of the periodic cor- 
rugation. In general, corrugated graphene is a perfect 



material with respect to geometrical effects due to the 
scale of the geometric interaction Vg ~ /m;f|k|, where 
for k w 100" 1 mn -1 (naturally formed ripples) equals 
Vg ~ 6 meV which is measurable. We believe that our 
results will lead to experimental verification of the phe- 
nomena suggested here. 
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